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We show that a special type of orbits, which are non- 
periodic and complex, 'Saddle Orbits' (SOs) describe accu- 
rately the quantal and experimental current oscillations in 
the Resonant Tunneling Diode in tilted fields. This is the first 
demonstration that one needs to abandon the periodic orbits 
(POs) paradigm in favour of this new and peculiar type of 
orbit. The SOs solve the puzzle of broad regions of experi- 
mental oscillations where we find no real or complex PO that 
can explain the data. The SOs succeeds in regimes involving 
several non-isolated POs, where PO formulas fail. 
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The Resonant Tunneling Diode (RTD) in tilted fields 
has recently been intensively investigated as an experi- 
mental probe of 'quantum chaos' It is widely con- 
sidered to be a paradigm of Periodic Orbit (PO) Theory 
in a real system, yet to date no PO formula has been 
shown to provide a full and quantitative description of 
the current. Several approaches were presented recently 
, expressing the tunneling current in terms of peri- 
odic orbits. We demonstrated previously jl(J that they 
could reproduce qualitatively the voltage range of experi- 
mental features like period doubling of the current. How- 
ever, they only yielded reasonable agreement for the am- 
plitudes of current oscillations in specific regimes, namely 
the stable (torus-quantization) region and its opposite ex- 
treme, the isolated unstable periodic orbit regions. They 
failed in an intermediate regime spanning a broad range 
of field values. There one finds regions where there is no 
real PO or alternatively competing non-isolated POs. As 
we argue below, it is not simply a question of improving 
the PO theory with uniform approximations. In the re- 
gions where there is no real PO, even complex 'ghosts' 
ijTTf POs cannot explain the experimental oscillations (for 
convenience we continue nevertheless to refer to these re- 
gions as 'ghost regions'). 

The main difficulty in establishing a semiclassical the- 
ory stems from the initial state (f>o describing the elec- 
trons prior to tunneling: one cannot simply exclude 
it from stationary phase considerations imposed on the 
rapidly varying function e lS ( z ' z ^ n of the classical action 
S of a trajectory. We show here that the correct station- 
ary phase condition arising from the theory proposed in 
H yields orbits of a new type - which we call 'Saddle Or- 
bits' (SOs). The SOs are non-periodic and complex, and 
describe the current accurately even in regimes where 
the previous semiclassical PO theories failed. They also 



show good agreement with the experimental amplitudes 
obtained from Bell Lab data The SOs are quite dis- 
tinct from the real closed orbits identified in atomic pho- 
toabsorption from localized ground states jl2). As for 
ghosts, the imaginary component of the action provides 
a damping term. We show that for SOs a further weight- 
ing term due to </>o can partially cancel this damping. 
This yields contributions which can decay slowly with 
decreasing h, also solving the puzzle of oscillations which 
persist far into the ghost regions. 

We recall briefly the RTD model JjJ . An electric field 
F (along x) and a magnetic field B in the x — z plane 
(at tilt angle 9 to the x axis) are applied to a double 
barrier quantum well of width L = 1200A. Electrons 
in a 2DEG accumulate at the first barrier and tunnel 
through both barriers giving rise to a tunneling current 
/. In the process they probe the classical dynamics - 
regular or chaotic- within the well. The current oscillates 
as a function of applied voltage V. After rescaling ^ 
with respect to B, the dynamics at given and ratio of 
injection energy to voltage (R — E/V ~ 0.15 for the 
Bell Lab experiments) depends only on the parameter 
e = V/LB 2 . 

The theoretical scaled current used for our quantal cal- 
culation and the starting point of the semiclassical theory 
is expressed as a density of states weighted by a tunnel- 
ing matrix element: 1(B) = Wj8(B — Bi). We used 
the Bardeen matrix element [J13|| form for Wi, which is 
an overlap between (f>o and the wave function ybi in the 
quantum well. In the experiments, incoherent processes 
such as phonon emission damp the current by e~ T l r , 
where T is the period of the motion of an electron in 
the well and r ~ 0.11 ps. Details of the experimental 
data reduction and quantum calculations were given in 
H[l0]]. We can obtain reliable experimental amplitudes 
in regimes where the current is dominated by a single 
frequency (pure period-one or period- two). 

For the semiclassics, one can re-express the Bardeen 
matrix element in terms of energy Green's functions and 
use their semiclassical expansion over classical paths to 
get the following expression for the tunneling current H : 

1(B) ocKe jdz [dz'Y,™i2 /2 z lS{z ' Z,) ~ BcOS(>{z2+Z ' 2)n (1) 

cl 

where m\2 — -J^— is an element of the monodromy ma- 

OPz 

trix niij, for a trajectory (x — 0,z) — > (x' — 0,z') 
which must connect both walls of the well. We consider 
here that the initial state is in the lowest Landau state: 
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4>o{z) = \J B cos 9/ir cxp(— B cos 9z 2 /2). The stationary 
phase condition applied to Eq.(Q) gives iff — B cos9z = 
= i-fip — BcosOz'. The contributing trajectories 
(z,p z ) (z',p' z ) will therefore satisfy the condition: 



Pz 



iB cos 9z 



Pz 



-iB cos 9z' 



(2) 



One sees clearly that these trajectories, which we call 
'Saddle Orbits' (SOs), will invariably be complex and 
non-periodic. It follows that we cannot consider the rep- 
etitions of a given SO -as one does for POs when inves- 
tigating period-doubled oscillations. Instead, one finds 
other SOs of quite different shapes with longer periods. 
We found that all the SOs which give a substantial con- 
tribution to the current have z = z' . They retrace them- 
selves once because of an intermediate bounce normal to 
a wall or a 'soft' bounce on the energy surface. 

In ll the B cos 9 term in Eq. (||) was neglected in order 
to get a formula in terms of POs with null momentum 
p z = = p' z . One can establish a link between a self- 
retracing SO and its PO counterpart by expressing the 
SO in term of an expansion around the PO. Making a 
Taylor expansion of the action and neglecting f^f for 
n > 2, one finds : 



zso 



zpo 



(3) 



where 5 = i cos 9 f ™^ 1 and fhij is the scaled monodromy 
matrix of the PO. We found that non self-retracing SOs, 
which correspond to segments of POs, are not relevant 
to these experiments. 

We plot some of these trajectories in Fig. [j]a)-d) , to- 
gether with their corresponding PO counterpart. The 
link between SOs and POs is evident in Fig. |l] a), which 
shows a regime where the SO and PO are very similar. 
In this case the x — z path of the SO is not very com- 
plex (the imaginary part is an order of magnitude smaller 
than the real part). The connection between SO and PO 
is generally not so apparent though. Fig. [j] b) shows the 
stable PO to, its related 'primitive' SO to-SO as well as 
another SO, which seems completely unrelated to them. 
This SO (2to-SO) plays, in fact, the role of the second 
repetition of the 'primitive' SO: its complex action and 
the real part of its period are twice (within 1%) those 
of to-SO. Despite the obvious difference in their path, 
some of the properties of the (SO-PO) pair such as their 
monodromy matrix or action can be very similar. 

Perhaps the most interesting situation is seen in Fig. |l| 
c), where the SO 'interpolates' between the ghost PO and 
the real PO. This fact is illustrated in Fig. [I] e), where we 
plotted the evolution of the starting position z (at x = 0) 
against e. One sees that a single SO is linked to a large 
number of POs. The latter are the basic traversing 2- 
bounce POs (to, t' , ..) to which have been attributed (in 
previous work) the period-one oscillations of the current. 
Their complicated dynamics has been well studied |D 



They undergo an infinite cascade of tangent bifurcations, 
where they disappear leaving a ghost. Subsequently at 
some lower e a similar PO reappears from the opposite 
edge of the Surface of Section. Over some e interval this 
re-entrant PO coexists with the ghost of the old PO. 
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FIG. 1. a)-d) Shape in x — z plane of the real part of 
SOs, with the related POs. a) 6 = 27°, e = 2000. Differences 
between the SO and its counterpart PO (Si) are minimal, b) 
8 — 11°, e = 20000. We show the main period-one PO to, its 
related SO as well as the SO which is related to its the second 
repetition 2t . c) 9 = 11°, e = 3000. The SO is between 
the real PO t' and the ghost PO t . d) 6 = 27°, e = 5000. 
A ghost PO (S") and a real PO {S") are present. The real 
part of the SO is closer in shape to the real part of the ghost, 
c) 9 — 11°. Evolution of the starting z (at x — 0) with e 
for the main period-one PO to, showing part of its infinite 
cascade of tangent bifurcations where to disappears, leaving 
a ghost while a similar PO appears at a lower e. Also shown 
is the behaviour of the related SO. We see that a single SO 
'interpolates' smoothly between the successive disappearing 
and re-entrant POs. This illustrates the striking simplicity 
and power of the SO approach in comparison with POs. 

In comparison, the corresponding SO behaves in a 
much simpler way, 'interpolating' between the successive 
POs. For example for e > 10000 the SO is related to to, 
which disappears in a tangent bifurcation at e = 6500. 
The SO remains close to the to ghost down to e = 3000, 
where it veers away towards the new orbit t' Q which ap- 
peared from the edge at e = 4300. One sees this fact 
clearly in Fig. [j] c), where the SO is between the to 
ghost and t' . The same happens at 9 = 27°, where the 
3-bounce PO S' disappears in a tangent bifurcation at 
e = 7700, leaving a ghost, while a similar PO S" appears 
at e = 5500. As seen in Fig. |l| d), the SO is very close 
the ghost of S', but it will soon approach the new PO at 
lower e. Note that SOs never disappear in bifurcations 
as they are non-periodic, but rather when they 'miss a 
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bounce' on the emitter wall because of the voltage drop. 

The semiclassical current is given by the straightfor- 
ward Gaussian integration resulting from the stationary 
phase approximation applied on (jl]). After normalizing 
to the amplitude at 9 = 0°, the current due to one SO is 
approximated by a simple analytical formula: 

pB(iS — cos 

1(B) = Re — ==== (4) 

■J — cos #772,12 + ?^2i / cos 9 + 2imn 

where fj, is a Maslov index, S and rhij are the scaled ac- 
tion and element of the classical monodromy matrix of 
the SO. One can also use the expansion in Eq.(||) to ap- 
proximate the SO current by the related PO. Expanding 
the action S of the SO around the PO up to second order, 
one finds: 

B (iS-cos 6z 2 )+ifnr/2 

1(B) ~ jfte (5) 
■J— cos0mi2 + m2i/ cos9 + 2imn 

where S and z are the scaled action and starting position 
of the PO. This is exactly the PO formula presented in ||] 
and tested in . So we see that the PO formula will give 
accurate results provided that the higher derivatives of 
S are small and that the expansion of the SO around the 
PO is justified. This latter point is the most relevant, as 
it is the prime reason why the PO formula fails in certain 
regions. 

We show in Fig. |^ a comparison between quantal, ex- 
perimental and semiclassical amplitudes for period one 
and two currents. The corresponding comparison with 
the PO formula was presented in Jio| . 

Fig. |2| a) shows the period one amplitudes at 9 = 11° 
in the unstable and 'ghost region', while Fig. |2| b) shows 
the stable torus regime. We found JIo| that for the sta- 
ble (e > 7000) and the chaotic (e < 2500) regions PO 
theory (using to and t' ) gave good agreement. These 
are regions where the SO and PO current are almost 
equal, according to Eq. (|5|). But the intermediate re- 
gion (2500 < e < 7000) was a puzzle: the PO formula 
including the ghost failed to account for the quantal and 
experimental oscillations, by a factor of 3. Fig. || a) 
shows that the SO completely solves this problem, giv- 
ing accurate results over the entire range from the torus 
regime to the chaotic regime, including the 'ghost region'. 

The main reason why Eq.(^) fails in that region was 
illustrated in Fig. [j] a): the SO cannot be approximated 
by the ghost PO, as it is also related to the new t' . In 
this case it is not because the third derivative of S is 
large. 

The same happens for the S' ghost at 27° [Fig. |d)]. 
The SO describes very well the broad plateau of quantal 
and experimental amplitudes, while POs failed. Once 
again, this is because the SO interpolates between the 
ghost (which appears at e = 7700) and the new PO S" 
(which appears at e = 5500), as illustrated in Fig. [I] d). 



Therefore, the failure of the PO theory in these regions is 
due to the sequence of tangent bifurcations, each followed 
by a new re-entrant PO, [as illustrated in Fig. |]f)], which 
rules out a simple connection between one SO and one 
PO. 
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FIG. 2. Quantal, experimental and semiclassical ampli- 
tudes. The SO labels indicate which PO they are related to. 
Period one at 8 — 11° in the unstable/ghost region (a) and 
in the stable (torus-quantization) regime (b). £o-SO describes 
very accurately the quantal current, even in the ghost region 
(2500 < € < 7000) where PO theories fail, (c) Period two 
at 6 = 11°. The SO formula improves over the PO formula, 
giving very accurately the quantal maximum (e ~ 15000). Ex- 
perimental amplitudes were not obtained for this case since 
there was a strong period one beat in this region, (d) Period 
two at 9 = 27°. As in (c), the peak at low e is well described 
by both the SO and the PO formulae (Si). S'-SO gives the 
contribution to the very broad plateau where no real PO is 
present (5000 < e < 8000). 2i -SO is responsible for the 
current for e > 17000, with no overlap with S'-SO. 

The SO also solves the intricate superposition of two 
non-isolated (in action and phase-space localization) POs 
ftPf : the second repetition of to (which yields quantized 
torus states) and S" for e ~ 13000, 9 = 27°. Indeed, there 
is no overlapping region for the related SOs, which now 
describe accurately the current, as shown in Fig. || d). 
In Fig. U c) one sees that the period doubling maximum 
around e ~ 15000 is described very precisely by the SO, 
an improvement over the PO results. Finally, we note 
that for both angles the period two maxima at low e 
~ 2000 are described equally well by either the PO or 
the SO formula. Here both trajectories (Si-SO and Si) 
are very similar as seen in Fig. [I] c). The experimental 
amplitudes are lower than theory, but this is consistent 
with a 10% uncertainty in r. 
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The strength and persistence of the contribution of 
these complex orbits, even in the region where the SO 
formula does not reduce to the PO formula (such as in 
the ghost regions) is quite remarkable. Usually (e.g., in 
the density of states or the photoabsorption spectra of 
atoms |ll],[n| ) , the contribution of complex ghost POs is 
extremely weak. They are exponentially damped away 
from the bifurcation (as e changes), since the imaginary 
part of the action increases as the ghost becomes more 
complex. However, in the RTD the tunneling amplitude 
includes an additional term due to the initial state: 



\I(B) 



oc e 



-B{Si+cos$(z%-zi)) 



(6) 



where the subscripts / and R denote the real and imag- 
inary parts. We see that the imaginary part Zi of the 
starting position of a SO can compensate for the damp- 
ing due to Si. This is the reason why even 'very' complex 
SOs can contribute. For instance, the S'-SO amplitude 
in the ghost region at 9 = 27° shows a broad plateau, 
and no exponential damping away from the bifurcation 
(e < 7700). 

Similarly, ghost POs are exponentially damped in the 
classical limit Ti — > (which corresponds in our scaled 
model to B — ► oo). In Fig. || we investigate the Ti — ► 
behaviour of the SO current in the two ghost regions. 
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FIG. 3. Semiclassical current 1(B) in two ghost regions 
showing the unusual Ti dependence, (a) 6 = 11°, e = 4000. 
The SO current is exponentially damped with B (i.e., in the 
classical limit), (b) 9 = 27°, e = 7000. The current is only 
linearly damped with B, showing a persistence in the classical 
limit never seen in ghost POs. 

In Fig. | a) (9 = 11°, e = 4000) the amplitude of t -SO 
is exponentially damped with B. This behaviour is seen 
in the experimental ghost region (see Fig. 5 of ||): along 
curves of constant V/B 2 the amplitudes decay rapidly 
with B. Fig. | b) (0 = 27°, e = 7000) shows a more 
surprising behaviour. Here the amplitude of S"-SO is not 
exponentially but linearly damped. This very surprising 
feature is seen experimentally (see Fig. 6 of ||). The 
explanation for the persistence of the complex SO in the 
classical limit is found again in 

We emphasize that this work is consistent with previ- 
ous studies of 'scarring' in the RTD. It has been found 



that quantum states localized near some isolated or mul- 
tiple POs can dominate the tunneling |||]]. We have also 
investigated wavefunctions and Wigncr distributions. We 
found that in the strong scarring regions, where the rele- 
vant scarring PO is only marginally unstable Q] (e.g. for 
9 = 27°, e = 10000), the real part of the SO is very close 
to the PO -within the quantum uncertainty. This can 
also be the case with scars carried by ghost POs (e.g. 
9 = 27° , e — 7000). This is reasonable since after all it is 
the bundle of classical trajectories in the neighbourhood 
of POs and SOs which scars or carries the electrons. Re- 
gions like the 9 = 11° ghost region where the SO and 
PO are really different in shape do not show strong scar- 
ring by single states and the quantal current is carried 
by broad clusters of states. 

We conclude that the SOs are a novel and success- 
ful way to approach semiclassical quantization of these 
types of chaotic systems. We recall that SOs arise solely 
from the inclusion of the initial state in the stationary 
phase condition. Hence one could expect SOs to be po- 
tentially relevant in the description of the expectation 
value of a quantal quantity (expressed as a density of 
states weighted by some matrix element Jl6| ) if the ob- 
servable in the matrix element is very localised and so 
varies as rapidly as e %s l n . 
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